Symmetry Energy, Temperature, Density and Isoscaling Parameter as a Function of 

Excitation energy in A ~ 100 mass region 
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The symmetry energy, temperature, density and isoscaling parameter, in ^^Ni + ^*Ni, ^®Fe + 
^*Ni and ®*Fe + ^*Fe reactions at beam energies of 30, 40 and 47 MeV/nucleon, are studied as a 
function of excitation energy of the multifragmenting source. It is shown that the decrease in the 
isoscaling parameter is related to the near flattening of the temperature in the caloric curve, and the 
decrease in the density and the symmetry energy with increasing excitation energy. The decrease 
in the symmetry energy is mainly a consequence of decreasing density with increasing excitation 
rather than the increasing temperature. The symmetry energy as a function of density obtained 
from the correlation is in close agreement with the form, Esym{p) = 31.6 (p/po)°'^^. 

PACS numbers: 25.70.Pq, 25.70.Mn, 26.50.+X 



Due to its vast implications ranging from how nucleons 
clusterize into nuclei at low densities to the structure and 
stability of neutron stars at high density, the interest in 
understanding the behavior of nuclear matter at temper- 
atures, densities and isospin (neutron-to-proton asymme- 
try) away from those of normal nuclear matter (T w 
MeV; po ~ 0.16 fm""^; N Z) has gained tremendous 
importance 0. 

Experimentally, the multifragmentation reaction 0, 0, 
^ IM1S> where a highly excited nucleus expands to a 
low density region and disassembles into many light and 
heavy fragments, provides the best possible means of 
studying nuclear matter at non-normal densities, tem- 
peratures and isospin. It has been shown from various 
experiments that the temperature as a function of ex- 
citation energy [caloric curve) in multifragmentation re- 
actions shows a near flattening, or a plateau-like region, 
at higher excitation energies. There are also indications 
that the density of the system decreases with increas- 
ing excitation energy Q . Recently, it has been shown 
that the isoscaling parameter, obtained from the frag- 
ment yield distribution, shows a decrease with increasing 
beam energy |lOlll]| . A decrease in the symmetry energy 
has also been experimentally observed fill ll^ . 

In this paper, we seek to understand the correlation be- 
tween the temperature, density and symmetry energy of 
a multifragmenting system as it evolves with excitation 
energy. Such a correlation is important for construct- 
ing the nuclear matter equation of state and studying 
the density dependence of the symmetry energy; a key 
unknown in the equation of state of asymmetric nuclear 
matter. 

We make use of the fragment yield distributions mea- 
sured mill El in s^Ni, sspe + ^^m, ^Spe reactions at 
30, 40 and 47 MeV/nucleon to determine the isoscaling 
parameter a, as a function of the excitation energy of the 
fragmenting source. The parameter a's were obtained 
by taking the ratio of the isotopic yields for two differ- 
ent pairs of reactions, ^^Fe -I- ^*Ni and ^*Ni -I- ^^Ni and 
58pg _^ 58pg j^^j 58]v^j _^ 58j^j described in Ref. 

The excitation energy of the source for each beam energy 



was determined by simulating the initial stage of the col- 
lision dynamics using the Boltzmann-Nordheim-Vlasov 
(BNV) model calculation 15]. The results were obtained 
at a time around 40 - 50 fm/c after the projectile had 
fused with the target nuclei and the quadrupole moment 
of the nucleon coordinates (used for identification of the 
deformation of the system) approached zero. These exci- 
tation energies were compared with those obtained from 
the systematic calorimetric measurements (see Ref. 0) 
for systems with mass {A 100), and similar to those 
studied in the present work, and are in good agreement. 
Fig. 1 shows the experimental isoscaling parameter a, 
as a function of the excitation energy obtained in the 
present study for both the pairs of reactions. A sys- 
tematic decrease in the absolute values of the isoscaling 
parameter with increasing excitation energy is observed 
for both pairs. The a parameters for the ^^Fe -|- ^^Fe and 
^*Ni -I- ^^Ni are about twice as large compared to those 
for the 5^Fe -I- ^^Ni and ^^Ni -t- ^*Ni pair of reactions. 
It is interesting to note that the difference in Z/A, i.e. 
/S.{Z/Af, of the composite systems in the ^*Fe -I- ^*Fe 
and ^^Ni -I- ^^Ni pair is also twice as large compared to 
the ^^Fe + ^^Ni and + pair. 

The experimental isoscaling parameter was compared 
to the predictions of the Statistical Multifragmentation 
Model (SMM) H |ll| calculations to study the depen- 
dence on the excitation energy and the isospin content. 
The initial parameters such as, the mass, charge and ex- 
citation energy of the fragmenting source for the calcu- 
lation, was obtained from the BNV calculations as dis- 
cussed above. To account for the possible uncertainties 
in the source parameters due to the loss of nucleons dur- 
ing pre-equilibrium emission, the calculations were also 
performed for smaller source sizes. The break-up density 
in the calculation was taken to be multiplicity-dependent 
and was varied from approximately 1/2 to 1/3 the sat- 
uration density. This was achieved by varying the free 
volume with the excitation energy as shown in Ref. Q. 
The form of the dependence was adopted from the work 
of Bondorf et ai, ^ (and shown by the solid curve 
in Fig. 4). It is known that the multiplicity-dependent 
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FIG. 1: Experimental isoscaling parameter a, as a function of 
excitation energy for the Fe + Fe and Ni + Ni pair of reaction 
(inverted triangles) , and Fe + Ni and Ni + Ni pair of reactions 
(solid circles) for the 30, 40 and 47 MeV/nucleon. The solid 
and the dotted curves are the statistical multifragmentation 
model calculations as discussed in the text. 
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FIG. 2: Temperature as a function of excitation energy for 
the Fe + Fe and Ni + Ni pair of reaction (inverted triangles), 
and Fe + Ni and Ni + Ni pair of reactions (solid circles) for 
the 30, 40 and 47 MeV/nucleon. The solid stars correspond 
to the measured values and are taken from Ref. 0. The solid 
and the dotted curve corresponds to the Fermi-gas relation. 



break-up density, which corresponds to a fixed interfrag- 
ment spacing and constant pressure at break-up, leads 
to a pronounced plateau in the caloric curve |5|, Il7j . A 
constant break-up density would lead to a steeper tem- 
perature versus excitation energy dependence. We will 
return to this point later in this paper. 

The symmetry energy in the calculation was varied un- 



til a reasonable agreement between the calculated and the 
measured a was obtained. It has been shown (Q.iiaLil9l|. 
that the symmetry energy in the statistical model calcu- 
lations is related to the isoscaling parameter through the 
relation, 

a^""' = ^[iZ/A)l-{Z/A)l] (1) 

where a^""^, is the isoscaling parameter for the hot 
primary fragments, i.e., before they sequentially decay 
into cold secondary fragments. Zi, Ai and Z2, A2 are 
the charge and the mass numbers of the fragmenting sys- 
tems. T is the temperature of the systems and Csym is 
the symmetry energy. In the above equation, the entropic 
contribution to the symmetry free energy is assumed to 
be small (the contribution becomes important at densi- 
ties below 0.008 fm~^ 20]), the symmetry energy can 
therefore be reliably substituted for the free energy. 

Fig. 1 shows the comparison between the SMM cal- 
culated and the measured a for both pairs of systems. 
The dotted curves correspond to the calculation for the 
primary fragments and the solid curves to the secondary 
fragments. The width in the curve is the measure of the 
uncertainty in the inputs to the SMM calculation. It 
is observed that, within the given uncertainties, the de- 
crease in the a values with increasing excitation energy 
and decreasing isospin difference A{Z/A)'^, of the sys- 
tems is well reproduced by the SMM calculation. One 
also notes that the effect of sequential decay effect on 
the isoscaling parameter is small as has been observed in 
several other studies using statistical models 0, 0| . 

We show in Fig. 2, the temperature as a function of 
excitation energy {caloric curve) obtained from the above 
SMM calculation that uses the excitation energy depen- 
dence of the break-up density to explain the observed 
isoscaling parameters. These are shown by the solid and 
inverted triangle symbols. Also shown in the figure are 
the experimentally measured caloric curve data compiled 
by Natowitz et al. 7] from various measurements for 
this mass range. The data from these measurements are 
shown collectively by solid star symbols and no distinc- 
tion is made among them. The Fermi-gas model predic- 
tions with inverse level density parameter Ko — 10 (solid 
and dashed lines), is also shown. It is evident from the 
figure that the temperatures obtained from the SMM cal- 
culations are in good agreement with the overall trend of 
the caloric curve. Somewhat lower value for the tempera- 
ture is observed when the break-up density of the system 
is kept constant at 1/3 the normal nuclear density. By al- 
lowing the break-up density to evolve with the excitation 
energy, a near plateau that agrees with the experimen- 
tally measured caloric curves is obtained. This assures 
that the input parameters used in the SMM calculation 
for comparing with the data are reasonable. 

The symmetry energies obtained from the statistical 
model comparison of the experimental isoscaling param- 
eter a, are as shown in Fig. 3. A steady decrease in the 
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FIG. 3: Symmetry energy as a function of excitation energy 
for the Fe + Fe and Ni + Ni pair of reactions (inverted trian- 
gles) , and Fe + Ni and Ni + Ni pair of reactions (solid circles) 
for the 30, 40 and 47 MeV/nucleon. 
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FIG. 4: Density as a function of excitation energy for the Fe 
+ Fe and Ni + Ni pair of reactions (inverted triangles), and 
Fe + Ni and Ni + Ni pair of reactions (solid circles) for the 
30, 40 and 47 MeV/nucleon. The solid stars correspond to 
those taken from Ref. Q. The open triangles are those from 
Ref. 0. The solid curve is from Ref. 13- 



symmetry energy with increasing excitation energy is ob- 
served for both pairs of systems. Such a decrease has also 
been observed in several other studies 0, 0, [S^, 113 . 
We have also estimated the effect of the symmetry en- 
ergy evolving durin g th e sequential de-excitation of the 
primary fragments [ll El. These are reflected in the 
large error bars shown in Fig. 3. 

The phase diagram of the multifragmenting system is 
two dimensional and hence the excitation energy depen- 



dence of the temperature (the caloric curve) must take 
into account the density dependence too. Often this de- 
pendence is neglected while studying the caloric curve. 
In the following, we attempt to extract the density of the 
fragmenting system as a function of excitation energy. It 
has been shown by Sobotka et al. that the plateau 
in the caloric curve could be a consequence of the ther- 
mal expansion of the system at higher excitation energy 
and decreasing density. By assuming that the decrease 
in the breakup density, as taken in the present statisti- 
cal multifragmentation calculation, can be approximated 
by the expanding Fermi gas model, and furthermore the 
temperature in Eq. 1 and the temperature in the Fermi- 
gas relation are related, one can extract the density as a 
function of excitation energy using the relation 



(2) 



In the above expression, the momentum and the fre- 
quency dependent factors in the effective mass ratio are 
taken to be one as expected at high excitation energies 
and low densities studied in this work [2^ l27l . 

The resulting densities for the two pairs of systems are 
shown in Fig. 4 by the solid circles and inverted trian- 
gles. For comparison, we also show the break-up densities 
obtained from the analysis of the apparent level density 
parameters required to fit the measure caloric curve by 
Natowitz et al. and those obtained by Viola et al. Q 
from the Coulomb barrier systematics that are required 
to fit the measured intermediate mass fragment kinetic 
energy spectra. One observes that the present results 
obtained by requiring to fit the measured isoscaling pa- 
rameters and the caloric curve are in good agreement 
with those obtained by Natowitz et al. The figure also 
shows the fixed freeze-out density of 1/3 (dashed line) 
and 1/6 (dotted line) the saturation density assumed in 
various statistical model comparisons. 

It is evident from figures 1, 2, 3 and 4 that the de- 
crease in the experimental isoscaling parameter a, sym- 
metry energy, break-up density, and the flattening of the 
temperature with increasing excitation energy are all cor- 
related. One can thus conclude that the expansion of the 
system during the multifragmentation process leads to 
a decrease in the isoscaling parameter, decrease in the 
symmetry energy and density, and the flattening of the 
temperature with excitation energy. Table I shows the 
correlated quantities obtained using Eqs. 1 and 2 for 
both pairs of systems. 

From the above correlation between the symmetry en- 
ergy as a function of excitation energy, and the density as 
a function of excitation energy, we obtain the symmetry 
energy as a function of density. This is shown in Fig. 5. 
for both pairs of systems. The temperature in the present 
work remains nearly constant for the range of excitation 
energies studied, the observed decrease in the symmetry 
energy with increasing excitation energy is therefore a 
consequence of decreasing density. This is also supported 
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TABLE I; Primary isoscaling parameter, temperature, density and symmetry energy obtained for Fe + Fe and Ni + Ni pair 
(top row), and Fc + Ni and Ni + Ni pair (bottom row) of systems at various excitation energies using Eqs. 1 and 2. 
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FIG. 5: Symmetry energy as a function of density for the Fe 
+ Fe and Ni + Ni pair of reaction (inverted triangles), and 
Fe + Ni and Ni + Ni pair of reactions (solid circles) for the 
30, 40 and 47 MeV/nucleon. The solid square corresponds 
to those from Ref. [33. The solid curve is the dependence 
obtained in Ref. |32L l3.?||. 

by microscopic calculations which shows an extremely 
slow evolution of the symmetry energy with temperature 
0, ll^. The evolution is practically negligible for the 
temperature range studied in this work. Also shown in 
Fig. 5 is the symmetry energy value of Khoa et al. [s^l , 
(solid square) obtained by fitting the experimental dif- 
ferential cross-section data in a charge exchange reaction 
using the isospin dependent optical potential. The solid 
curve corresponds to the dependence, ^sym{p) = 31.6 
(p/po)^'^^, obtained by comparing the present data with 
the Antisymmetrized Molecular Dynamic (AMD) calcu- 
lation [31| , in previous work [ifl, Is^, 123 • The similarities 
in the density dependence of the symmetry energy ob- 
tained from the present statistical model approach and 
the AMD model approach is intriguing. It should be 
noted that the symmetry energy shown by the solid curve 
in the figure relates to the volume part of the symmetry 



energy as in infinite nuclear matter, whereas, the symme- 
try energy obtained from the present study (solid circles 
and inverted triangles) relates to the fragment that are fi- 
nite and has surface contribution. The similarity between 
the two can be understood in terms of the weakening of 
the surface symmetry free energy when the fragments are 
being formed. During the density fluctuation in uniform 
low density matter, the fragments are not completely iso- 
lated and continue to interact with each other, resulting 
in a decrease in the surface contribution as predicted by 
Ono et al. M^. The present observation therefore lends 
credence to the fact that it is possible to directly obtain 
the properties of infinite nuclear matter from the frag- 
ments produced in the multifragmcntation process 3^. 
More studies are required to illustrate this point further. 

In summary, we have studied the isoscaling parameter, 
symmetry energy, temperature and density as a function 
of excitation energy in reactions populating A ~ 100 mass 
region. It is observed that the decrease in the experimen- 
tal isoscaling parameter a, symmetry energy, breakup 
density, and the fiattening of the temperature with in- 
creasing excitation energy are related to each other. The 
observed decrease in the symmetry energy with increas- 
ing excitation energy appears to be mainly a consequence 
of decreasing density. The symmetry energy as a function 
of density obtained from the present study is consistent 
with those obtained from the dynamical AMD calcula- 
tion, indicating that the surface contribution to the sym- 
metry energy could be small. 
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